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Abstract 



We include the relativistic lattice KP hierarcliy, introduced by Gibbons and Kupersh- 
midt, into the r-matrix framework. An r-matrix account of the nonrelativistic lattice 

C^ , KP hierarchy is also provided for the reader's convenience. All relativistic construc- 

tions are regular one-parameter perturbations of the nonrelativistic ones. We derive 
in a simple way the linear Hamiltonian structure of the relativistic lattice KP, and 

f~^ ' find for the first time its quadratic Hamiltonian structure. Amasingly, the latter turns 

\^ , out to coincide with its nonrelativistic counterpart (a phenomenon, known previously 

^D ' only for the simplest case of the relativistic Toda lattice) . 

o\ 

1 Introduction 

One of the basic objects in the theory of integrable lattices (or differential-difference 
_^ systems) is the so called lattice KP hierarchy Q, |6[. This is a system of an infinite number 

c/3 . of commuting flows on the phase space consisting of an infinite number of fields 



j^ ■ The simplest ("first") flow of this hierarchy is governed by the following equations of 



motion: 






(i+i) Aj+i)^ ^"^'^^ 



A proper language for such infinite-dimensional systems is the differential-difference cal- 
culus developed in |Q, ^. However, all essential properties of this hierarchy hold also for 
its finite-dimensional reductions described as follows. 

First of all, we shall consider here only systems with a finite number of fields. We obtain 
an ?TT,-field system, if we set a^ = for all j > m. We shall use the notation TL^ for this 
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reduced system, because it is a direct generalization of the Toda lattice, which appears in 
this framework as the first nontrivial case TL = TL2 with two fields. 

Further, we restrict ourselves here to the case of periodic or open-end boundary con- 
ditions, when all fields contain only a finite number A^ of variables, so that in (^]^) one 
should replace A; € Z through 1 < k < N. More precisely, in case of periodic boundary 
conditions all subscripts in (|1.2|) are supposed to belong to Z/iVZ, and in case of open-end 
boundary conditions we enforce 

6fc = for k < 1 and for k > N, 

(i) 
a^ =0 for k <1 and for k > N — j. 

One of the outstanding features of this hierarchy is integrability; another one is its bi- 
Hamiltonian structure. What we want to address in the present paper, is an r-matrix 
theory as the (probably most direct) route towards understanding these both properties. 
Actually, the r-matrix interpretation of the lattice KP hierarchy is a more or less 
established thing nowadays ||l|, |ll|, 15, p. What is really new in this paper, is a similar 



interpretation of the relativistic lattice KP hierarchy. This appeared i n pO| , Q as a natural 
generalization of the relativistic Toda lattice, invented by Ruijsenaars |l2|. The relativistic 
ansatz of Q (referred to later on as "the first construction") leads to a one-parameter 
perturbation of the lattice KP hierarchy. An interesting "splitting" phenomenon takes 
place: to each (polynomial) flow of the lattice KP hierarchy there correspond two flows of 
its relativistic conterpart, one of them being still polynomial, and another one rational in 
coordinates. It turns out that the rational relativistic perturbation of the "first" flow ( |1.2[) 
remains nice and elegant: 

4" «£■ 



1 + abk+i 1 + abk-i' 

(1.3) 

(i) _ Ji) ( ^k+j bk \ , ( a,k "fc-i \ 



"'^ "'^ V 1 + abk+j l + abkj yi + abk+j+i 1 + abk-i J ' 

while the polynomial relativistic perturbation is described by a system of messy equations, 
right-hand side of each one depending on all fields. In the m-field situation we shall 
denote the flow (|1.3| ) and the whole hierarchy attached to it by RTLm (a). It is a direct 
generalization of the "minus first" flow of the relativistic Toda hierarchy. 

It is known that this relativistic lattice KP hierarchy consists of commuting flows, and 
also one Hamiltonian structure was known for it previously 0, ^, ^]. In the present work 
we find also the second compatible Hamiltonian structure, and, moreover, give an r-matrix 
account for both. Needless to say, that this automatically implies also the integrability of 
this hierarchy. Amasingly, while the first Hamiltonian structure of the relativistic hierar- 
chy is a perturbation of its nonrelativistic counterpart, the second Hamiltonian structures 
of the both hierarchies literally coincide! 

The relativistic ansatz of Q leads to Lax equations of a somewhat non-standard type. 
We introduce a certain (invertible) gauge transformation, such that the transformed hier- 
archy possesses the standard Lax representations. In other words, it may be considered 
as included into the lattice KP hierarchy. It is a highly non-trivial situation, if one takes 
into account that the lattice KP may be seen also as a limiting case of the relativistic 
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lattice KP! This gauge transformed relativistic hierarchy (referred to as "the second con- 
struction") also contains two perturbations, a polynomial one and a rational one, for each 
nonrelativistic flow, but this time the polynomial flows are simpler. For example, the 
polynomial "first" flow in the second relativistic construction reads: 



bk = {l + abk) (4') - a^^\ 



n. _ n. I 6fe+j - 6fc + aflfc+j - ««fc-l j + \^k - "fc-1 ) ■ 



(1.4) 



We denote this flow, as well as the whole corresponding hierarchy, in the m-field situation 
as RTLm (ce); it serves as a direct generalization of the "first" flow of the relativistic Toda 
hierarchy. 

Remarkably, in the case m = 2, i.e. of the relativistic Toda, both constructions lead to 
one and the same hierarchy, but this is no more the case for m > 2\ 

Also the second relativistic lattice KP turns out to be bi-Hamiltonian and integrable; 
we give in this paper an r-matrix interpretation of one of its Hamiltonian structures. 

The paper is built as follows. In Sect, g-g we recall several notions and results from 
the r-matrix theory (those of Sect. ^ being to a certain extent new). Sect. ^ is devoted 
to fixing the basic notations. The lattice KP hierarchy and its linear and quadratic 
invariant Poisson structures are dealt with in Sect. 0-|^. The similar discussion of the 
first relativistic construction is contained in Sect. 11-1^ (so that the presentation of really 



new results starts in Sect. |12| ). Further, we introduce the second relativistic construction 
in Sect. 14, ^, and discuss its Hamiltonian properties in Sect. |T^. As illustrations, we 



provide concrete results for the three-field systems in Sect. |10|, |17|. A final discussion takes 



place in Sect. 18. 



2 Linear r-matrix Poisson structure 

We first recall the (by now classical) construction of the linear r-matrix structure on the 
dual to a Lie algebra. Let g he a Lie algebra, carrying a nondegenerate scalar product 
(•, ■), allowing to identify the dual space g* with g. In what follows we shall not distinguish 
between g and g*, but try to use the letters X, Y, Z consequently for the elements of g, 
and the letters L, M for the elements of g* . Suppose that the scalar product is invariant 
with respect to the Lie bracket in g, i.e. 

{X,[Y,Z]) = {[X,Y],Z) yx,Y,zeg. 

Recall that the gradient V(/7 : g >-^ g of a smooth function ip on g is defined via 

y M eg. 



de 



e=0 

Recall also that Ad-invariant functions <p on g are characterized by the property 

[Vip{L),L]=0 VLSff. (2.1) 

(Of course, we mean here actually Ad*-invariant functions on g* , characterized by 
ad*V^(L) • L = 0; 

we shall consequently avoid such remarks from now on). 
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Definition 2.1 [13] Let Yv he a linear operator on g. A linear r -matrix bracket on g 



corresponding to the operator R is defined by: 

{<^,^}i(L) = ^([R(Vv9(L)),VV^(L)] + [Vv9(L),R(VV^(L))],L). (2.2) 

If this is indeed a Poisson bracket, it will he denoted by PBi(R). 



Theorem 2.2 |13| A sufficient condition for (2.<:) to define a Poisson bracket is given 



by the modified Yang Baxter equation for the operator R; 

[R{X),R{Y)]-R{[R{X),Y] + [X,R{Y)]) = -a[X,Y] yX,Y e g, (2.3) 

where a is an arbitrary constant. This equation is denoted mYB(R;a). 



Actually, ( p.3| ) is a sufficient condition for 

[X,Y]n = l{[R{X),Y] + [XMy)]) (2.4) 



to define a new Lie bracket on g, and then (2.2) is nothing but a Lie-Poisson bracket 



on g* corresponding to this new Lie algebra structure on g. 

One of the most important properties of the linear r-matrix bracket is the following one. 



Theorem 2.3 |13| a) Hamiltonian equations of motion on g corresponding to an Ad- 



invariant Hamilton function ip, have the Lax form 

L=[L,C], C = ^R{V^iL)). (2.5) 

b) Ad-invariant functions on g are in involution with respect to the bracket PBi(R). 

3 Generalized linear r-matrix Poisson structure 

Let now g have an additional structure of an associative algebra, with the multiplication 
{X, Y) I— > X ■ Y. The standard definition of the Lie bracket on g is then [X, Y] = 
X -Y — Y ■ X. Let g carry a bi-invariant scalar product, the property expressed by the 
equalities 

{X, Y ■Z) = {X ■Y,Z) = {Z-X,Y) yX,Y,Z£g 

(this assures the invariance of the scalar product in the previous sense relevant to the Lie 
algebra structure). We again identify g* with g with the help of this scalar product. 

Actually, there exist infinitely many ways to turn g into a Lie algebra, not only the 
standard one mentioned above. Indeed, fix an arbitrary element J^ £ g, and define the 
corresponding Lie bracket as 

[X,Y]^^^ = X ■J'-Y -Y -J^-X. (3.1) 

We can replace the Lie bracket [•, •] in all constructions of the previous Section by the Lie 
bracket [•, ■]^'^', and come in this way to the following definition. 
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Definition 3.1 Let R be a linear operator on g. A generalized linear r-matrix brac- 
ket on g corresponding to the operator R and the element J- is defined by: 

{<^,^}i(L) = i([R(Vv.(L)), VV^(L)](^) + [Vv9(L),R(VV^(L))](-^),L). (3.2) 

If this is indeed a Poisson bracket, it will be denoted by PBi(R;.7-'). 



Of course, a sufficient condition for ( |3.2[ ) to be a Poisson bracket is that the following 
expression defines a new structure of a Lie algebra on g: 

[x,y]if) = \ ([R(x),y](^) + [x,R(y)](^)) . (3.3) 

A sufficient condition for this is given, clearly, by the mYB(R;a), where the Lie bracket 
[•, •] is replaced through [•, •]^'^\ i.e. 

[R(x), R(y)](^)-R ([R(x), y](^) + [X, R(y)](-^)) = -« [x, y](-^) v x, y g £,.(3.4) 

Notice that Ad-invariant functions ip of the Lie algebra {g, [•, •]' ') are characterized by 
the equation 

L-Vi){L)-T-T-V^{L)-L = {) yL£g. (3.5) 

The analog of Theorem |2.3| can now be formulated. 

Theorem 3.2 a) Let ip be an Ad-invariant function of the Lie algebra (gi, [•, •]' ') . Then 
the Hamiltonian equation on g generated by the Hamilton function ip in the Poisson bracket 
PBi(R;J^) reads: 

L = L-C2-Ci-L, (3.6) 

where 

Cx = \t- R(VV^(L)), C2 = ^R(VV'(L)) • T. (3.7) 

b) Ad-invariant functions of the algebra {g,[-,-] ) o,re in involution with respect to 
the bracket PBi(R;J^). 

The above notions simplify under some additional structural assumptions. Let g be an 
associative algebra with unit, and let T ^ g be an invertible element. Then: 

Proposition 3.3 a) The function t/j is an Ad-invariant of the algebra (g,[-,-]^-^^), i.e. 
( |3.5| ) is fulfilled, if and only if 

ij{L) = (^ (L • .F-i) = ^ (^-1 • L) , 

where (p is an Ad-invariant of the algebra g with the standard Lie bracket. In this case 

V^{L) = T-^ ■Vip{L- T-^) = Vip {T-^ • L) ■ J^-\ 

b) The modified Yang-Baxter equation (3.4) in the Lie algebra [g, [■, •]' ') for the oper- 



ator R is equivalent to the modified Yang-Baxter equation ^2.3i ) in the algebra g with the 
standard Lie bracket for either of the two linear operators on g, 

Yli{X)=r-R{r-^ -X) and R2(X) = R (X • J^^^) • J^. (3.8) 

c) Under the map L i— > L-T^^ the bracket PBi(R;^) is pushed to PBi(Ri), and under 
the map L i— >• T~^ ■ L the bracket PBi(R;jF) is pushed to PBi(R2). 
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4 Quadratic r-matrix Poisson structure 

Let again g be an associative algebra with a nondegenerate invariant scalar product. Let 
Ai, A2, S be three linear operators on g, Ai and A2 being skew-symmetric: 

Al = -Ai, A^ = -A2. (4.1) 



Definition 4.1 [15] A quadratic r-matrix bracket on g corresponding to the triple 



Al, A2, S is defined by: 

W,i^}2{L) = \{Ki{d'^{L)),d'^{L)) - \{K2{d^{L)),dij{L)) 

(4.2) 
+^(S(ML)),dV(^)) - ^(S*(dV(L)),d^(L)), 

where we denote for brevity 

dip{L)=L-Vip{L), d'ip{L) = \/ip{L)-L. (4.3) 

If this expression indeed defines a Poisson bracket, we shall denote it by PB2(Ai, A2, S). 

In what follows we usually suppose the following condition to be satisfied: 

Al + S = A2 + S* = R. (4.4) 

Then a linearization of PB2(Ai,A2,S) in the unit element of g coincides with PBi(R), 
and we call the former a quadratization of the latter. 



Theorem 4.2 |15| A sufficient condition for (4-^) to be a Poisson bracket is given by the 



equations (1^4-A ) ^''^^ 

mYB(R;a), mYB(Ai;a), mYB(A2;a). (4.5) 

Under these conditions the bracket PB2(Ai,A2,S) is compatible with PBi(R). 

If the operator R is skew-symmetric and satisfies mYB(R; a), then the Poisson bracket 



PB2(R, R, 0) is called Sklyanin bracket [13|. The brackets PB2(A, A, S) with a skew- 



symmetric operator A and a symmetric operator S were introduced in [^, |l^] . 

One of the most important properties of the r-matrix brackets is the following one. 



Theorem 4.3 |15| Let the condition (4^) be satisfied. Then: 

a) Hamiltonian equations of motion on g corresponding to an Ad-invariant Hamilton 
function ip, have the Lax form 

L = [L,C], C = hi{d^{L)). (4.6) 

b) Ad-invariant functions on g are in involution with respect to PB2(Ai, A2, S). 
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5 Quadratic r-matrix structures on direct products 

Quadratic r-matrix brackets have interesting and important features when considered on 

n 

a "big" algebra g = (H) 9- This algebra carries a (nondegenerate, bi-invariant) scalar 
product 

n 

{{L,M)) = Y,{Lk,Mk). 

k=l 

Working with linear operators on g, we use the following natural notations. Let A : g i— s- g 
be a linear operator, let {A(L))- be the ith component of A{L); then we set 

n 

(A(L)), = ^(A),,(L,). (5.1) 

i=i 

For a smooth function ^{L) on g we also denote by Vj<I>, dj^, d',^ the jth. components 
of the corresponding objects. We define also the monodromy maps Tj : g >-^ g , 1 < j < n, 
by the formula 

Tj{L) = Lj ■ . . . ■ Li ■ Ln ■ ■ ■ ■ ■ Lj+i. (5.2) 

Now let Ai, A2, S be linear operators on g satisfying conditions analogous to (fl.lj) 
and to ( [4.5[ ). One has, obviously: 

((Ai),,)* = -(Ai),„ ((A2).,)* = -(A2),., ((S).,r = (S*),.. 
Then one can define the bracket PB2(Ai, A2, S) on g. In components it reads: 

{<^,^}2{L) = -Y, ((Ai),,(4<I>),d^vi/) _ - ^ ((A2).,(d,$),d,M/) 

i,j=l i,j=l 

(5-3) 



+i ^((S).,(d,<I>),d^*) -i 5^((S*).,(4ci>),d,,I.). 



2 ^^ '''"'"''' ' 2 



Theorem 5.1 |16| Let g he equipped with the Poisson bracket PB2(Ai, A2, S). Suppose 
that the following relations hold: 

(Ai),+ij+i + (S),+ij = {A2)j,j + (S*)jj+i = R for all 1 < j < n; (5.4) 

(Ai)i+ij+i = -(S)j+ij = (S*)j,j+i = -(A2)ij for i / j. (5.5) 



Then each map Tj : g h^ g fjS.^ ) is Poisson, if the target space g is equipped with the 
Poisson bracket 

PB2((Ai),+ij+i,(A2),j,(S),+ij). 

Hamilton function of the form ^{L) = ip{Ln • . . . • Li), where (p is an Ad-invariant function 
on g, generates Hamiltonian equations of motion on g having the form of Lax triads: 

Lj = LjCj^i - CjLj, Cj = -R{d<p{Tj)). (5.6) 

(In all formulas the subscripts should be taken (mod n)). 



418 Yu.B. Suris 

This theorem is a far-reaching generahzation of the corresponding result for the Sklyanin 
bracket PB2(R, R, 0) on g, which arises when S = 0, Ai = A2 = R, and (R)y = R6ij. In 
this case each map Tj : g ^^ g (5^) is Poisson, if the target space g is equipped with the 
Sklyanin bracket PB2(R, R, 0) [13[. Certain generalizations of the latter result appeared 
also, e.g., in [^, ^, |l^, but in all previously known formulations only few nonvanishing 
"operator entries" for the operators Ai, A2, S were allowed, namely "diagonal" ones for 



Ai, A2, and "subdiagonal" ones for S. In other words, all operators in (5.5) had to vanish. 
In the present paper, however, we shall need this Theorem in its full generality (albeit only 
for n = 2). For the first time this general bracket was applied in [^]. 

We have discussed above the r-matrix origin of Lax equations, such as (|2.5| ), ( |4.6D , 
or (^]^) . If one is concerned with a Lax representation of a certain Hamiltonian flow 

x = {R,x} (5.7) 

on a Poisson manifold {X, {■,■}), then finding an r-matrix interpretation for it consists 
of finding an r-matrix bracket on g (or on g) such that the Lax matrix map L : X >-^ g 
(resp. L : ^ 1-^ g) is a Poisson map. Then the manifold consisting of the Lax matrices is 
a Poisson submanifold. 

6 Notations 

6.1 Algebras 

Two concrete algebras play the basic role in our presentation. They are well suited to 
describe various lattice systems with the so called open-end and periodic boundary con- 
ditions, respectively. Here are the relevant definitions. 

For the open-end case we always set g = gl{N), the algebra oi N x N matrices with 
the usual matrix product, the Lie bracket [X, Y] = XY — YX, and the nondegenerate 
bi-invariant scalar product {X,Y) = tr{XY). As a linear space, g may be represented as 
a direct sum 

N-l 

9= gp, 

p=-N+l 

where the subspace gp consists of matrices 
X = 2_^ XjkEjk- 

j-k=p 

(Here and below Eji^. stands for the matrix whose only nonzero entry is on the intersection 
of the jth row and the kth. column and is equal to 1). The following sets are subalgebras 
oi g: 

N-l 

p=i 

so that, for instance, g^Q consists of lower triangular matrices, g^Q consists of strictly 
upper triangular matrices, and g^Q consists of diagonal matrices. We shall always set 

9+ = g>0: 9^=9<o- 



N-l 

®9p, 


9<o -- 


-1 
- 9p, 


9<o -- 



= 9p 


p=0 




p=-N+l 




p=_Af+l 
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Notice that, as a linear space, 

9 = 9+®9 , 

In the periodic case we always choose gf as a certain twisted loop algebra over gl{N). 
A loop algebra over gl{N) is an algebra of Laurent polynomials with coefficients from 
gl{N) and a natural commutator [XA-', yA^] = [X, yjA-'^^. Our twisted algebra g is a 
subalgebra singled out by the additional condition 

g = {X(A) G gl{N) [A, X-^] : nX{X)n-^ = X{ioX)} , 

where Q, = diag [l,uj, ■ ■ ■ jCO^"^^, uj = exp(27ri/A^). In other words, elements of g satisfy 

^w = E^' E ^t^^^- (6.1) 

(mod AT) 

The nondegenerate bi-invariant scalar product is chosen as 

(X(A),y(A))=tr(X(A)y(A))o, (6.2) 

the subscript denoting the free term of the formal Laurent series. 
As a linear space, g is again a direct sum 

oo 

9= 5p, 

p=— oo 



where gp consists of matrices 
X = X^ 2_^ XjkEjk- 



j — k^p 
(mod AT) 



We have the subalgebras 

9>o = 9p, 9>o = 9p, 9<o = 9p, 9<o = £/ 

p>0 p>0 p<0 p<0 

Again, we set 

9+=9>o, 9-=9<o, 

so that, as a linear space, 

9 = 9+®9-- 



p- 
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6.2 Operators 

In both the open-end and the periodic cases we shah denote by P>o, P>Oi etc. the pro- 
jections from g onto the corresponding subspace g^Q, g^Q, etc., along the complementing 
subspace. For two such projections we use special notations: 

7r+ = P>o, 7r_ = P<o. 

The basic operator governing the hierarchies of Lax equations, is 

R = 7r+-7r_. (6.3) 

An important property of this operator is given by the following general statement: 

Theorem 6.1 |j3, |13| Let g as a linear space be a direct sum of its two subalgebras: 

9 = 9+® 9-, 

and let ■7t± denote projections from g onto g^ along the complementary subspace. Then 
the operator 

R = 7r+ — 7r_ 

satisfies the modified Yang-Baxter equation mYB(R, a) with a = 1. 

Let us remark that for the operators ( |6.3D the Lax equations ( |2.5| ), ( |4.6| ) are equivalent to 

L=[L,B] = [A,L], where B = 7r+{f{L)), A = Tr_{f{L)) (6.4) 

with f{L) = Vip{L) and f{L) = dip{L), respectively. Similarly, the Lax equations on the 
"big" algebra g (^) for the operator (|6.3D are equivalent to: 



^j = ^jBj-i - BjLj = AjLj - LjAj_i, (6.5) 

where 

B,=7T+{fiTj}), Aj=7T^{f{T,)). 

We shall call such Lax equations standard ones. 

Denote by Rq, Pq the skew-symmetric and the symmetric parts of the operator ( |6.3D , 
respectively: 

Ro = (R-R*)/2, Po = (R + R*)/2. 

It is easy to see that 

Ro = P>o - P<o, 

and the notation Pq agrees with the previously defined projection to giQ- ^^ the open-end 
case Pq assigns to each matrix X its diagonal part, and in the periodic case Pq assigns to 
each Laurent series X{X) its free term. 

Let the skew-symmetric operator W act on gfg according to 

N 
j<k j>k j=l 
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where 




Wkj = sgn{k - j) = < 0, k = j (6.7) 



and extend W on the rest of g according to W = W o Pg. Finally, define: 

Ai = Ro + W, Aa = Ro - W, S = Pq - W, S* = Pq + W. (6.8) 

These operators will be basic building blocks in all quadratic r-matrix brackets appearing 
in this paper. 

7 Multi-field analog of the Toda lattice 

The {m + l)-field phase space is, in the periodic case, 

T„+i=M(™+i)^(6,a«,...,a(-)), (7.1) 

where the vectors 

b={h,...,bN) and a(^) = (4^\...,ai^)) (i = l,...,m) 



represent the m + 1 fields. In the open-end case the fields a'-^' consist oi N — j variables 
(^a['\...,a^^l.) only. 

Consider the Lax matrix 

N N m N 

T (b, a^'\. . . , a^"^), a) = a j; Ek+i,k + ^ bkEkk + J^ A'^^ J^ a^^^Ek^k+j- (7.2) 

fc=l fc=l j=l k=l 

The subspace of g consisting of all such matrices will be denoted 

m 

T^+i=f e05_j-. (7.3) 

N 

(We introduced here the notation £ = X ^ E^+i^k', recall that in the periodic case all 

k=l 

subscripts are considered (mod N), so that Ejy^i^^ = Ei^j\[). 

In the open-end case the Lax matrix, by convention, is given by 

N-l N m N-j 

T (b, a(i), . . . , a(")) = Y, Ek+i,k + Y. b^Ekk + Y.Y. ^k^Ek,k+j, 

k=l k=l j=l k=l 

and the space of all such matrices is denoted still by ( |7.3| ). (So, in the open-end case 

Af-l 
£ = Yl Ek+i,k)- In what follows we shall formulate the results and give the proofs 
k=l 

mainly for the periodic case, since for the open-end case they appear to be parallel and 
only more simple. 
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Proposition 7.1 [Q] The Lax equation 

t=[T,B], (7.4) 

with the Lax matrix ( \7.30 and the auxiliary matrix 

N N 

B (b, a(i) , . . . , a^'") ,X)=Y1 ^kEkk + A ^ E^+i^u, (7.5) 

fc=i fc=i 

is equivalent to the following system of differential equations: 

4'^ = ^f (h+, - bk) + (4'^'^ - «£'0 ' l<j<m-l, (7.6) 

V '^fc ~ ^k \"k+m — Ok) ■ 

Proof - an elementary matrix calculation. I 



The system (|7.q ) is what will be called TL^+i - the generalized (m + l)-field Toda 
lattice. It serves as a direct generalization of the usual Toda lattice TL. A consistent 
notation for TL would be TL2, i.e. it corresponds to the value m = 1. 

8 Linear r-matrix structure for TL^+i 

Obviously, the matrix B from (|7.5| ) is nothing but the projection 

B = 7r+{T) (8.1) 

in the standard decomposition of the algebra g. So, the Lax representation of the system 
TLm+i is of the standard form ( |6.4D with f{T) = T. We now give an r-matrix interpre- 
tation to this Lax representation, defining simultaneosly the whole hierarchy TL^+i. We 
start with the linear r-matrix Poisson structure on g. 



Theorem 8.1 |13, 11 1 The set T^+i is a Poisson submanifold in the algebra g equipped 



with the linear bracket PBi(R). 

Proof. We have to demonstrate that at each point T G T^+i every Hamiltonian vec- 
tor field {^p,-}i is tangent to T^+i. To do this, notice that the value of the general 
Hamiltonian vector field in the point T is given by: 

2{v9,r}i = [T,R{Vip)]+R*{[T,V^]). 
Obviously, this can be represented in two equivalent forms: 

{<^,r}i = [T,P>o(Vv.)] -P>o([r,Vv9]), (8.2) 

and 

{if, r}i = -[T, P<o(Vv9)] + P<o([T, W^]). (8.3) 
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1 
Now take into account that T G ® 9j- Then, according to ( |8.2[) , the value of {(p,T}i 

belongs to Qj^ and, according to ( ^.3D , this value belongs also to g^. Hence it 

i>—m j<0 



belongs to Qj, which is the tangent space to the manifold T^+i in g. I 

j=-m 

It remains to calculate the induced bracket on 7^+i. In what follows we always define 
Poisson brackets by writing down all nonvanishing brackets among the coordinate func- 
tions. 

Theorem 8.2 The coordinate representation of the restriction o/PBi(R) to the suhma- 
nifold Tm+i is given by the formulas: 

{6.,4^)}^ = -4-'), {a?,6.+,}^ = -a?, (8.4) 

{4^4i}^ = -«r'^ (^ + J<H. (8.5) 

Proof. We have, obviously: 

Vbk = R(V6fc) = Ekk, V4*^ = R (v4'0 = >^'Ek+i,k. (8.6) 



Hence 

\bk, «£ I = {T, A-' [Ekk, Ei-f-j/]) = (T, X^ E£^j^i{6k/+j — Sk/)) = a^ [^k,t+j — ^k/), 
which coincides with ( |8.4| ). Analogously, 

l^k '"£ I = {'^y'^^^'^[Ek+i,k,Ei+j^i]) = {T,X^^^Ek:+i/5k/+j - Ei+j^k^l^k+i) 

— a^ Ok/+j — a J. Oi^k+i, 

which coincides with (^.5|). One shows similarly that the brackets {bk,bj}i vanish. This 
finishes the proof. I 

Corollary. ^ Th 

(8.5), with the Hamilton function 

N N 



Corollary. [|4[ The system TL^+i ( |7.q j is Hamiltonian with respect to the bracket ( \8.4) , 
with the Hamilton function 

N N 



fc=i fe=i 



Theorem |8 . 1| not only provides us with the Hamiltonian structure for the system TL^+i, 
but gives also an r-matrix explanation to the Lax representation from Theorem 7.1, and 
allows to define the TL^+i hierarchy as a set of flows 

T=[r,vr+(Vv9(r))] 

with Ad-invariant functions ip. 

Also the complete integrability of TL^+i with respect to the linear Poisson bracket can 
be studied on the base of Theorem 3.1, because it guarantees involutivity of all spectral 
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invariants of the matrix T (b, a'^^ , • • • , a*-™"^ , A) . Therefore, the complete integrabihty will 
follow as soon as it is proved that the number of independent spectral invariants is large 
enough. This is almost obvious in the case m = 1, i.e. for the usual TL, but for m > 1 
it is no longer evident, and, moreover, is no longer true in general. However, it may be 
shown that in the periodic case one gets the sufficient number of independent spectral 
invariants. In the open-end case the situation is somewhat more delicate. One has to 
introduce integrals which are not spectral invariants |^, ^. 

9 Quadratic r-matrix structure for TL^+i 



Theorem 9.1 ||15| , ^ The set T^+i is a Poisson submanifold in the algebra g equipped 
with the quadratic bracket PB2(Ai, A2, S). 

Proof. The Hamiltonian vector field on g with an arbitrary Hamilton function 99 with 
respect to the bracket PB2(Ai,A2,S) is given by: 

2{ip, T]2=T- AiiVifT) - A2(rV(/7) -T + T- S{TVip) - S*{VipT) • T. (9.1) 

For the operators ( |6.8D this can be represented in two equivalent forms: 

2{ip,T}2 = 2rP>o(V(^T) - 2P>o(TV^)r 

+rPo(Vv9r) - Po(TV99)r + TPo{TVip) - Po(Vv9r)r (9.2) 

+rw(Vv9r) + w{TVip)T - rw(rv(^) - w(Vv9r)r, 

and 

2{^,T}2 = -2rP<o(V99r) + 2P<o(TVvp)T 

-TPoiVipT) + Po(rV(^)r + rPo(rV(^) - Po(v<^T)r (9.3) 

+rw(Vv9r) + w{TVip)T - rw(rv(^) - w(Vv9r)r. 

The first expression assures that for T E T^+i this vector field takes the value in ^ gr •, 

j>-m 

and the second expression yields that that this vector field belongs also to ^ g,:- Hence 

1 
it belongs to © gj. It remains to prove that the gi component of this vector field 

j=-m 

vanishes. Up to now everything held true for an arbitrary operator W with the values in 



gfQ. Now the concrete expression (6.6) for W becomes crucial. From ( [7.2[ ), ( |9.3| ) we have 

(9.4) 



the following expression for the component in question: 

2Pi ({v9, r}2) = -SPoC^^T) + Po(rVv9)^ + £Po{TVip) - Po{VcpT)£ 
+£W{VipT) + W(rVv9)^ - £W{TVip) - W{VipT)£. 
Due to W = W o Pq the condition for ( |9.4D to vanish may be presented as 

W{D)£ - £W{D) + D£ + £D = W{D')£ - £W{D') + D'£ + £D', (9.5) 

where 

D = PoiTVif), D' = Po(V(^T). (9.6) 
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Now a direct calculation shows that for the operator W in ( |6.6[) and for an arbitrary 
diagonal matrix D one has: 

{2A tr(L>) El TV, periodic case 
(9.7) 
0, open — end case 



This proves ( |9.5D for the matrices (|9.6D , since ti{D) = tr{T'Vip) = tr(V(/9T) = ti(D') (note 
that in the open-end case the last argumentation is superfluos). I 



Theorem 9.1 gives the second Hamiltonian interpretation of the system TL^+i, and 
yields simultaneously the complete integrability of TL^+i in this Hamiltonian formulation 
(at least in the periodic case). Moreover, it delivers an alternative r-matrix explaination 
of the Lax representation from Theorem |7.1| . It remains to calculate the coordinate rep- 
resentation of the restriction of PB2(Ai, A2, S) to Tm+i- In the sequel we often use in our 
formulas the following notational convention: 

af^ = bk; 4?^ =0 for i < or i > m, (9.8) 

whenever applicable. 

Theorem 9.2 The bracket induced on T^+i by PB2(Ai, A2,S) is given by the following 
formulas (in ^9.13^ ) we suppose that i < j, and in ( 9.1^ , l \9.l3i) we suppose that i < j): 



{bk,bk+i}2 = -a'}^\ 




(9.9) 


{h,a^&.}, = -at'\ 




(9.10) 


{6.,a?}^ = -6.a?, 


{^^^\bk+,}^ = -a^^\+j, 


(9.11) 


{«i^^agm},=-«r^"^^ 


(9.12) 


\^k '"-fc /2 ~ "-k ^k ' 


{^\4l,-l = -^'4i,-. 


(9.13) 




-a^-'^a^;P il<i<^), 


(9.14) 




-ot^'o^;? U-^ + l<i<J)■ 


(9.15) 



Proof. Notice that without the convention ( |9.8D the second term on the right-hand side of 
( 9.14| ) would have to be multiplied by xe (li niin(i, m — j)), and, likewise, the second term 
on the right-hand side of ( |9.15| ) would have to be multiplied by xe (i ~ ^ + 1, ™iii(j, m — i)), 
where 

^^^"'^) = \0, else 

is the characteristic function of the interval [a,/3]. Notice also that under the conven- 
tion ( p^ ) the formulas (p^), (9.10) might be seen as particular cases of ( |9.12| ). Similarly, 
( 9.11| ) might be seen as a particular case of ( |9.13|) . 
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Toward the calculation of the induced bracket, we use the defining formula 

2{^,V'}2 = 2(p>o(V(/pr),vV'r) -2(p>o(rV(^),rvV') 
+(Po(V(/^T), WT) - (Po(rv^),rvv^) + (Po(tv^), vv^t) - (Po(V(^r),rvv) 
+(w(Vv9r), VV'T) + {w{TVip),TVyj) - (w(rvv9), VV'T) - (w(V(^r),rvV'). 

(9.16) 
Take here 

(^(r) = a«, V(T) = 4-5„ 

so that 
Further, we find: 

m 

V^T = \'+^Ek+^,k-l + Y. A*"''af £^fc+i,fc+/3, (9.17) 

m 

TVif = A*+ii?fe+,+i,fe + Y, A*"^afi-/j^fc+i-/3,A:. (9.18) 

We consider first the contribution to the Poisson bracket {ip, ^p}2 = \ aj^ , a^^i \ from the 
first line in ( 9.16| ): 



V+ii?,+,,fc_i + ^A*-^afSfc+.,fc+/j, Y. ^'^^''kleEk 



13=0 7=i+l 

k+e+j-y-^k+e+j-j,k+e 



-{X'+^Ek+i+l,k + Y^'' '^H+i-l3^k+i-P,k, Y ^^ ^"'k+l+j—t^k 
\ /3=0 7=J+1 

Calculating these scalar products, we find: 

m 

" Y \^k+l^-^Mj^i/+l - O'k+l+j--,^0/+3-l^i+'^r 
7=j+l 

i— 1 m 
/3=0 7=i+l 






where 



^iy^ = X<?(-i> ^ - J - l)xK^ -m,i-j-l)- xeil, m - j)xeil, i) 
= Xi (- min(j, m-i),i- j -1) -xe (1, min(i, m - j)) . 



(9.19) 
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Assuming, for the sake of definiteness, that i < j, we see that the intervals of the two 
characteristic functions in the last line do not intersect. So, we found the contributions to 
the Poisson bracket < cij^ ,aj^^£ \ described by the formulas (|9.9D, ( 9.10| ), ( |9.12| ), and the 
second terms on the right-hand sides of ( 9.14 ), ( 9.15| ). Now we turn to the contribution 



of the remaining part of (9.16). To this end, we notice that 



.W 



Po(V99r) = a):> Ek+^,k+^, Po(rV(^) = a):>Ekk. 



W; 



Using also ( |6.6D , we find that the contribution under consideration is equal to 






where 



'ki 



+Wk+i,k+e+j + Wk^k+e — Wk^k+e+j — Wk+i,k+i)- 



(9.20) 



(9.21) 



A direct analysis allows to find a compact expression for this coefficient. In the case i < j 
we have: 



'-ke 



Xe{-j,i-j) -Xi{0,i). 



(9.22) 



Again, by i < j the intervals of these two characteristic functions do not intersect, and 
hy i = j their only common point is ^ = 0, so that the last formula exactly describes the 
contribution into the Poisson bracket made by ( p. 11 ), ( p. 13 ), and the first terms on the 
right-hand sides of (9.14), ( 9.15| ). This finishes the proof. I 

Corollary. Q The system TL^+i ( 17. 6| ^ is Hamiltonian with respect to the bracket given 



in Theorem 9.2, with the Hamilton function 



N 



Hi 



Y,bk = {tiiT)),. 



(9.23) 



k=l 



10 Example: TL3, 

the three-field analog of the Toda lattice 

In order to illustrate the above results, we give their specialization for the case next in 
complexity after the usual Toda lattice TL, i.e. for m = 2. The equations of motion of 
the system TL3 read: 



a 



(1) „(i) 



fc-i' 



a(1) 



.(2) 



,(1) 



{bk+i -bk)+ (a 



,(2) .(2) 



k-1 ' 



(10.1) 



af {bk 



'Jk+2 -h) ■ 
The linear invariant Poisson bracket of this system is given by: 

{".,4"}, = -4". {4".'>»H-.}, = -" 



{hJ^>}, 



"k ' 



\cik'\h+2 






i 



-a 



k ' 

(2) 
k ' 



{ 






,(2) 



(10.2) 
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The quadratic invariant Poisson bracket, compatible with the previous one, is given by: 


{bk,bk+i}2 = -%\ [h 


'4^^}2 = -bka^k\ {ak\bk+i}^ = -a^^bk+i, 


{&fc,4+i}2 = -«?^ {° 


l^\&..2}^ = -4^\ {bk,a^^'}^ = -bka^^\ 




r . (10.3) 

\^k ''^fc+1/2 ~ fc "fc+i' 


f„(l) „(2) 1 _ _„(l)a(2) 


r„(2) ^(1) 1 _ _a(2)„(i) 

\"fc '"fc+2/2 ~ °fc "'k+2' 


/„(2) J2) 1 _ (2) (2) 


r„(2) „(2) 1 _ _.(2)„(2) 



11 Multi-field analog of the relativistic 
Toda lattice: the first construction 

The essence of the relativistic ansatz of Q consists in multiplying the nonrelativistic Lax 
matrix by J^''^{X), where 

N 

T{X) = I-aXj2Ek+i,k, (11.1) 

fc=i 

and the (small) parameter a is the relativistic parameter (the inverse speed of light). So 
the Lax matrix of the relativistic analog of TL^+i is 

TJ^-^ or J^-^T, 

where T is the matrix from (|7.2| ). It will be convenint to consider also the matrices 

Ti (6, o(i), . . . , a("^). A) =I + aT {b,a(^\. . . ,a(™), A) • .^-^(A) 

= £(6,o«,...,a('"),A)-.^-i(A), 

T2 (b, a(i), . . . , a^™), a) = / + aT-\X) ■ T (b, a'^^\ . . . , a^'^l x) 

^ ^ ^ ^ (11.3) 

= ^-i(A)->C(6,aW,...,a("^),A), 

where 

£ {b, a(i), . . . , a^"^). A) = T{X) + aT [b, a^^\ . . . , a^™). A) 

^ "" ^ r^ (11-4) 

= J](l + abk)Ekk + « XI ^"^ IZ 4^^fc,fc+i- 
fc=i j=l fe=i 

The arising hierarchy, as well as its simplest ("minus first") flow will be denoted 
RTL)^^-^(a). We shall use the notation 

7^Tj,;\ = m(-+i)^ (6,a(i), . . . .a^™)) (11.5) 

for the phase variables of this hierarchy. 
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Proposition 11.1 |^ Consider the matrix differential equation 

t = TB2- BiT, (11.6) 

with the auxiliary matrices 



N , N 

k=l " ' " fc=l 



E— — JT-E^fcfc + A^ — — j-E^fc+i,fc, (11-7) 

1 + abu ^-^ 1 + abu 

^ b ^ 1 
'S2 = V '^Ekk + XY^ ; Ek+ik- (11-8) 

fc=i '^ k=i "^^ 

The matrices Bi satisfy the relation 

TB2 - BiT = 0, (11.9) 

which assures that the equation ( |11.6D is equivalent to either of the following two Lax 
equations: 

fi = [Ti,Bi], i = l,2. (11.10) 

These matrix differential equations are equivalent to the following system: 

^k "'k-l 



1 + abk+1 1 + abk-i ' 

^(i) ^ ^(i) f ^k+j ^J^_\ j^ I ^k^ "fc-i I ^ (11-11) 

^ ^ \l + abk+j 1 + abkJ yl + abk+j+i 1 + abk-i J ' 

1 + abk+m 1 + ah 



<' = <' ( -r-r^ -r-r^ ) > l < i < m - l. 



In 1^] this system was called "the shadow relativistic lattice KP" . Obviously, the system 
(11.11) is a one-parameter perturbation of ( [7.6[ ). In order to define the hierarchy attached 



to the flow RTL^^^(a) and to find its Hamiltonian structure, we shall use the r-matrix 
theory. 

12 Linear r-matrix structure for RTL^^j^(a) 

Let us consider the generalized r-matrix bracket PBi(R;.7^) introduced in Sect. ^, with 
the standard operator 

R = 7r+-7r_. (12.1) 

First of all, we have to be sure that this is indeed a Poisson bracket. Notice that the 



matrix T is invertible. Therefore the constructions of Proposition 3^ are applicable. 
Lemma 12.1 Let the linear operators Ri^2 on g be defined by the formulas 

Ri{X)=T ■R{T~^X) , R2{X)=Yl{XT-^)-T. (12.2) 

Then these operators satisfy the modified Yang-Baxter equation. 
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Proof. The proofs for Ri and R2 are similar, therefore we restrict ourselves to the Ri 
case. We prove that this operator can be represented as 

Ri = P(+)-P(_), (12.3) 

where P(±) denote the projections from g to the corresponding subspaces g(±) in a certain 
decomposition of g (as a linear space) into a direct sum 

a = a(+)®a{-). (12.4) 

each one of the subspaces gij^\ being also a subalgebra and a Lie subalgebra of g. Then 
the statement of the Lemma is assured by the general construction of Theorem 6J.. To 
prove the above claim, set 

and 

P(+)(X) = ^7r+ {J^-^X) , P(_)(X) = J^TT^ {T-'X) . 



Obviously, from these definitions there follows the direct decomposition ( 12.4 ), the repre- 
sentation ( |12.3| ), as well as the following statements: 

P(+) (X) + P(_) (X) = X, P(+) (X) G »(+) , P(_) (X) G 5(_) . 

It remains to show that gi^\ are subalgebras. But this follows from the obvious relations 

XTYGg+ yX,YGg^ ^ TX ■ TY e g(^+) V .FX, .FY € »(+) , 
and 

XTY£g_ yX,Y£g^ =^ TX ■ TY € g(^_) V J^X, .FF G f/(_) . 

The Lemma is proved. I 

Remark. It is important to notice that the operators Ri and R2 correspond to different 
splittings ( |12.4|) . In both cases one has gf_^_\ = g_^ = g>Q, but gi\ are different subalge- 
bras of gi<o for the two cases. Let us give also somewhat more explicit formulas for the 
operators Ri,2- We use the following notation for the negative part of X: 

N 

7r_(X) = J]A-^ j;xi-')i?fc,fc+,. (12.5) 

l>0 fc=l 

Then for the operators Ri^2 = P(+) ~ '^i-) ^^ have: 

P(+)(X) = 7r+(X)+ai,2(X), P(_)(X)=7r_(X)-ai,2(X), (12.6) 

respectively, where the operators o"i^2 ■ 9 ^^ 9o ^ct according to the formula 

N N 

ai(X) = ^a'^xi-'^Ek+,,,+e, ct2{X) = Y,a'Y,^^k'^ E^k- (12.7) 

i>0 fc=l i>0 k=l 

In words: the operators cri^2 assign the weights a^ to the elements of g_i, and shift all 
entries of 7r_(X) to the diagonal positions: cri along the rows where they stand, and (T2 
along the columns where they stand. In this form the operators o"i^2 appeared in pH, R] 
(in the latter reference - under the name "relativistic symbol"). 
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Theorem 12.2 a) The set T^+i is a Poisson submanifold in g equipped with the linear 
r -matrix bracket PBi(R;^). 

b) Let if be an Ad-invariant function on g, and let il^{T) = ip (TT^^^ = tp (^J^^^T^. 
Then the Hamiltonian equation on g with respect to PBi(R, J^) with the Hamilton function 
ip{T) reads: 

t = TC2- CiT, (12.8) 

where 

Ci = ^Ri (V^ {TT'')) = ^-^ • R {J'-'^V {TT-')) , (12.9) 

C2 = ^R2 (Vv? {T~'T)) = Ir {Vif {T-'T) J^-') ■ T; (12.10) 

this equation may be properly restricted to T^+i- The evolution of the matrices 7i = 
/ + aTJ^^^ and T2 = I -\- aJ^~^T is described by the usual Lax equations 

ti = [Ti,Ci], f2 = [T2,C2]. (12.11) 

c) For two Ad-invariant functions (^1, ip2 on g the functions ipiiT) = ipi (TJ^^^^ = 
ifi (^~^T) and ip2{T) = '^2 {TJ-^^^ = ip2 [J-~^T^ are in involution with respect to 
PBi(R,.F). 



Proof follows closely the arguments used in Theorem 8.1. We denote by {•,•}! the 
Poisson bracket PBi(R; J^) on g. For an arbitrary Hamilton function ■0 on g, the value of 
the Hamiltonian vector field {'0, -ji in the point T G g is given by: 

2{^, T}i = TR(VV')J^ - J^R(VV')T + R* {TVipT - TV^T) . 

This may be represented in the following two alternative forms: 

{V, r}i = rp>o(vv).F - .FP>o(vv^)r - p>o {tv^t - tv^jT) , (12.12) 

and 

{V, r}i = -rp<o(vv^).F + .FP<o(vv)T + p<o (rvv.^ - j^v^t) . (12.13) 

1 

Now take into account that T G fl'j ^^'^ ^ ^ 9o®9i- Then, according to ( |12.12| ), the 

j=-m 

value of {'il^,T}i belongs to ^ gj. From the first glance, from ( 12.13 ) one can derive 

j>-m 
only that {tlj,T}i belongs to ^ gj, but, rewriting ( 12.13|) as 

i<i 

aii/j, r}i = -/:p<o(vV').f + .^p<o(vv)/: + p<o (cv^jJ^ - j^v^jC) , 

we see that actually the stronger inclusion {V', T}i G ^ g^ holds. Hence, the value of the 

i<o 



vector field under consideration belongs to gr,-, which means that Tm+i is a Poisson 

j=-m 

submanifold in g. Other statements of the Theorem follow from the general constructions 
of Sect. H ■ 
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In the formulation and in the proof of the next Theorem we use the conventions (|9.8| ) . 

Theorem 12.3 The bracket induced on T^+i by PBi(R, J^) has the following coordinate 
representation: 

{bk,bk+i}i = aa^k\ (12.14) 

{6.,ag,}^ = a4^+^), {4^\6.+,+i}^ = «4^-+^ (12.16) 

{4^4i}, = -«r\ (12.17) 

Proof. To calculate the induced bracket, we denote 

^{T) = a^\ V'(T)=ag„ (12.19) 

so that 

V99 = RiVif) = X'Ek+i,k, VV' = R(VV') = A^'^fc+^+i,fc+£. (12.20) 

We find the following value for the Poisson bracket {^, ip}i = \ o^ , dj^^i \ '■ 



T, A* ^ Ek+i,k • il — aX 2_^ En+i^nj ■ E^+i+j^k+e 

N 

T, A* ^ Ek+t+j,k+e. ■ (l — a\ 2_^ En+i,n] ■ Ek+i,k 

n=l 

T, y^^ [di^-jEk+i^k+i — h,iEk+i.+j,k) 



—aX" ^ (5£ _j_i-E'fc+j^fc+£ — ^i,i+iEk+e.+j,, 

This proves the Theorem. I 

The Theorems above deliver a Hamiltonian interpretation of the flow ( 11. ll]) , as well as 
an explanation of the Lax representations of this flow given in Proposition [11. 1| . Indeed, 
as a Corollary be find the following statement. 
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Proposition 12.4 ^, ^ The flow ( 11.11 ) on TZT^^_^^ is Hamiltonian with respect to the 
linear bracket { 12. 1^ - ^12. 18[) , with the Hamilton function 



N 



Hn 



-a 



'^log(l + a6fc). 



(12.21) 



fc=i 



Of course, this function is singular in a, but this can be repaired by adding the following 
Casimir function: 



N 



N 



a-iHi = Q-iJ]6fc + ^ 



(1). 

k ' 



(12.22) 



fc=l k=l 

obviously, the resulting Hamilton function is regular in a, and, moreover, has an asymtotics 

N N 

2 



^ N N 

Ho + a-^Hi = j; 6| + ^ 4^) + 0{a). 



k=l 



k=l 



The integrals of motion (12.21), (12.22) may be represented as the following functions 
of T: 

Ho = ^ {T-^T) = 99 {TT-^) = -a-2 (log det(Ti,2))o , Hi = (tr(ri,2))o • (12.23) 

For the Hamilton function Ho in ( |12.23| ) we find: 



This explains the expressions ( p,1.7 ), ( 11.^ ) for the auxiliary matrices B\^2'- indeed, up to 
adding a matrix proportional to /, 



i3i 



-oT^T -Ti^^L-^), 



Bo 



-a-^TT+ (C-^) ■ T, 



in agreement with (|2j), (pIlo| ). The formulas ([l2^ , ( |T2^) , ( pTTo|) define various 
members of the hierarchy RTL,^j^]^(a). The equations of motion of all other members, 
except ( |11.11| ), are much more long and complicated. The reader should verify that this 
is the case already for the Hamilton function iIj{T) = tp (^"^T) such that S/ip{T) = T, 
which corresponds to the "first" flow of the hierarchy, as opposed to the "minus first" 
flow RTL^^]^(q:). As an example, we give here the formulas of this "first" flow for m = 2: 

b, = il + ab,)(af^ ^« .^^^(^) ^.^(^) 



'k - Ofc-i + ""fc - "«fc-2 ; ' 



■M) 



M) 



,(2) 



,(1) 



(1) 



.2^(2) 



bk+1 -bk + aOfc+i - a«fc-i + " "fc+i + " 
(2) 



0^(2) 0^(2) 



a a 



2.(2) 



fc-1 



a a 



k-2 



+a):> (1 + abk+i) - a^"_!i (1 + abk) , 



.(2) _ (2) /, , ,(1) (I) ,^2„(2) ,^2„(2) 0(2) ^2^(2) \ 

Ofc - H ["k+2 - Ok + aof^^^ - aaf^_-^ + a a^^^ + " H+i " " H-i " " H-2) ■ 
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13 Quadratic r-matrix structure for RTL^^^(q:) 

Now we prove a very remarkable Theorem. According to it, the quadratic invariant Poisson 
bracket of the usual Toda lattice is at the same time invariant with respect to the flows 
of the relativistic Toda hierarchy. 

Theorem 13.1 Supply g with the quadratic r-matrix bracket PB2(Ai,A2,S) defined by 
the operators (6.8). Let f be an Ad-invariant function on g, and consider the Haniiltonian 
flow on g with the Hamilton function 

ipiT) = LP (TT~^) = ^ [T-^T) . (13.1) 

The equations of motion of this flow read: 

t = TC2- CiT, (13.2) 

where 

Ci = ^Ri {dip (TT-^)) = ^-^ • R [J'-^d^ [TT-^)) , (13.3) 

C2 = ^R2 {d^ {T-^T)) = Ir {d^ {T-^T) ^-1) • T. (13.4) 

Proof. For an arbitrary Hamilton function ip{T), the Hamiltonian equations of motion 
have the form ( 13. 2| ) with the matrices 

2Ci = A2 (TVV'(r)) + S* {Vi){T)T) , (13.5) 

2C2 = Ai (vV'(r)r) + s (Tv^(r)) . (13.6) 

For the functions ip{T) = tp (TT~^) with Ad-invariant (p we have, obviously: 

Vip{T) = T-^ ■ V<p [TT-^) = Vip {T~^T) ■ T'^ , 
so that 

TVip{T) = dip (TT-^) , Vtp{T)T = dip {J^-^T) . 

Hence, denoting 

X = T-^ ■ dip {TT-^) = dip (T'^T) ■ T-^, 
we can rewrite (|13.5| ), ( |13.6|) as 

2Ci = k2{TX) + S*(XJ^), (13.7) 

2C2 = kx{XT) + ^{TX). (13.8) 



We now prove the following two identities. For an arbitrary matrix X S gf let 

N 

,(-1) : 



P_i(X) = X_i = A-1 Y. x{ ^^Ek,k+i. (13.9) 
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Then 

MTX) + S*{XT) = TR{X) - 2axt^^/, (13.10) 

Ai(XJ^) + S{TX) = R{X)T - 2axS^^^/. (13.11) 

This will, obviously, imply the statement of the Theorem. Since ( 13.10D , (13.11) are proved 
analogously, we do this only for the first one of them. According to (|6.8| ), and J- = I — a£, 
we have: 

A2{TX) + S*{XT) = Rq{TX) + Po(^.^) - W(J-X) + W(X^) 

= R(X) - a (Ro(^X) + Vq{X£) - W{£X) + W{X£)) . 

The term proportional to a is equal to: 

P>o(f X) - P<o(^X) + Po(X^) + W ([X_i, f ]) 

= fP>o(X)-fP<_i(X)+X_i^ + W([X_i,^]) 
= £K{X) + £X^i + X_i^ + W {[X^i,£]) . 
Now a direct calculation shows that for an arbitrary matrix X_i € g_i there holds: 

f X_i + X^i£ + W ([X_i, ^]) = 2xS^^^/, 

and the proof is finished. (Notice that in the open-end case the left-hand side of the last 
equation simply vanishes). I 



The Theorem just proved yields that the matrix differential equations ( 13.2 ) can 
be properly restricted to an arbitrary Poisson submanifold of the quadratic bracket 
PB2(Ai, A2, S), and the resulting equations are Hamiltonian with respect to the restriction 
of this bracket. In particular, this holds for the Poisson submanifold T^+i, and we find 
the following statement. 

Theorem 13.2 The hierarchy RTL^_j_-^(a) is Hamiltonian with respect to the quadratic 
invariant Poisson bracket (9.{i)- (9AB^ of the non-relativistic hierarchy TL^+i. This brac- 
ket is compatible with the linear bracket of Theorem 12.S^ . 

In particular: 



Proposition 13.3 The system ( 11.11 ) is Hamiltonian with respect to the bracket {•,-}2 
on TZT^j^^ given by the formulas \9.(\ )-{ 9.15 ). The corresponding Hamilton function is 
equal to 



N 



-aHo = a ^ y^log(l + abk). 



k=l 
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14 Introducing the gauge transformed hierarchy 

The Lax equations of RTL^_j_j^(a) are of a somewhat nonstandard type. We now find 
a gauge transformation bringing the Lax equations into the standard form. To this end, 
we start with the matrices J^{X), C{a,b,X) for the RTL^^j^(— a) hierarchy, i.e. we change 
a to —a in the formulas ( p.l.lD , ( [11.41 ). As it fohows from Theorems 12.2] , 13.1, we have 
the following Lax representations for the flows of the RTL^^^(— a) hierarchy: 

fi = [Ti,Bi] = [A,Ti], i = l,2, (14.1) 

where 

Bi=T-7T+{T-'f{CJ^-')), Ai=T-7T^{T''f{£J'-')), (14.2) 

B2 = vr+ (/ {T-'C) T-^) -T, A2 = vr_ (/ {J^-^C) T'^) ■ T, (14.3) 

with an Ad-covariant function f ■ g >-^ g (we have f{T) = Vip{{T — I) /a) and f{T) = 
dip{{T—I)/a) for the Hamiltonian flows with the Hamilton function ip(TJ^~^^ = (/?(^^^T) 
in the linear and the quadratic bracket, respectively). 
Now introduce two arbitrary diagonal matrices 

^i = diag{uji,uj2,. ■ ■ ,ujn), il.2 = diag{uj[,uJ2, ■ ■ ■ ,uj']y). (14.4) 

It will be convenient for us to denote 

Qi ■ Ti ■ nr^ = Tr\ i = 1,2. (14.5) 

The gauge transformed Lax equations for the matrices Tj are easily derived from ( |14.1[) 
and ([T45|) : 

ti = [Ti,Bi] = [Ai,Ti], i = l,2, (14.6) 

where 

Bi = QiBiQ-^ - nin-\ Ai = QiAiQ-^ + tliQ-^ (14.7) 

^Prom the expressions ( |14.2D , (14.3) one sees immediately that 

Bi G g>Q, Ai e gf<o, 
and the same holds for the matrices Bi, Ai: 

B^ G £/>o, Ai G £/<o. (14.8) 

Our aim is to assure that 

B^€g>o = g+, Ai£g^o = g_. (14.9) 

To this end it is necessary only to kill the diagonal entries of Ai, i.e. to define the entries 
of the matrices Oj in such a way that 

ni^-^ = -Po(A). (14.10) 

So, all we need to do, is to determine the diagonal parts of the matrices Ai- 
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Lemma 14.1 // 

Po(^i) = diag(ai, . . .,3.^), P0M2) = diag(ai, . . .,a'^), 

then 

where C is some constant (independent on k). 
Proof. Consider the matrix equations 

^i-F - J^.42 = 0, C = AiC-CA2. (14.11) 

Due to Ai G g<Q and the definition of !F, we find from the first equation in ( 14.111) : 



a'fc = afe+i. (14.12) 

On the other hand, since C G g^Q, we find, considering the Qq part of the second equation 
in (114.111) : 



-a6fc = (1 -a6fc)(afc -a'fc) =^ a^+i - a^ 



abk 



1 - abk 

This implies the statement of the Lemma. I 

So, let us define the entries of the diagonal matrices Oj, i = 1,2, by the following 
formulas: 

fc-l k 

u^k = S"-' n(l - «^^)' ^fe = '^'"' 11(1 - ab,), (14.13) 

where we set in the periodic case 

N 

5'' = l[{l-ab,r\ 
i=i 

and in the open-end case simply 6 = 1. Of course, in the periodic case 5 is a spectral 
invariant of the Lax matrix, and therefore an integral of motion of an arbitrary flow of 
the hierarchy. 

Theorem 14.2 Fix an Ad-covariant function f ■ g ^^ g and consider the corresponding 
flow of the hierarchy RTL^^-^(— a). Then the evolution of the matrices Ti is described by 
the standard Lax equations 



fi = [r„vr+ (/ {Tr^))] = [vr. (/ {Tr^)) ,T,] , i = 1,2. 
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Proof. From Lemma 14.1 and the definitions ( |14.13| ) we find: 



afc — C , a^ — G y, 

which means that 

Po{Ai)+ninr^ = C-I. (14.14) 

Therefore, renaming Ai — C ■ I hy Ai, and Bi -\- C ■ I hy Bi (which does not affect the Lax 
equations ( |14.6| )), we come to the desired property ( 14.91 ). Further, from ( 14.2[ ), ( |14.3D we 
derive: 

Bi + Ai = f{%) 

and it follows that 

B, + A = f{Tr^). (14.15) 

Therefore Bi and Ai are with necessity the 7r+ and 7r_ projections of f(T~ ), respec- 
tively. ■ 

In particular, the flow RTL^_j_-^(— a) is characterized, as we have seen, by f{T) = 
a~^T^^. Therefore its gauge transformed Lax representation reads: 



T, = a-'[Ti,7r+{T,)]. 



(-) 



In other words, the image of the flow RTL„^_j_]^(— a) in the gauge transformed hierarchy 
becomes the "first" flow. 



15 Multi-field analog of the relativistic 
Toda lattice: the second construction 

An easy calculation shows that (we still consider £, T attached to the RTL^^-|^(— a) 
hierarchy) : 

fti ■ j^{x) -n^^ = lU c(^) , . . . , c^'") , 6X) , 

where 

N N 

L [d, c(i) , . . . , c^'") , a) = ^(1 + adk)Ekk + aXj2 Ek+i,k, (15.1) 

fc=i fc=i 

m N 

U {d, c(i), . . . , c('"), a) =I-aY,^''Y. ''k^^k,k+j, (15.2) 

i=i fc=i 
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and the variables d, c^^' are defined by the formulas 



B(a) : 4 



hk 



Si) 



Aj) 



l-abk' 



0(1- abk+i) 

i=0 



(15.3) 



We want now to study the gauge transformed hierarchy in its own rights. The phase 
space of this second version of the generalized relativistic Toda lattice with m + 1 fields, 
abbreviated RTL^_j_^(a), is, in the periodic case, the space 



7^r^+\ = M(-+l)^fd,c«,...,c('-) 



(15.4) 



Here 



d={di,...,dN) and c(-') = (ci^\...,cS^)) (j = l,...,m) 

are the m + 1 fields. The Lax matrix map (L, U^^) : TZT^^J^ i—fg®gis defined by the 
formulas ( |15.1D , ( |15.2| ). We set also 



Ti = LU-\ 



U-'L. 



(15.5) 



Proposition 15.1 The Lax triads 

L = LB2 - BiL, U = UB2 - BiU 



Bi (d,c^'\. . . ,c('"\ a) =Y,(dk + acl^l,) Ekk + X^Ek+i,k, 



with the Lax matrices (15.1), (15.i) and the auxiliary matrices 

N N 

.(1) 
-A;-] 

fc=l k=l 

N N 

B2 {d,c^'\. . . ,c('"\ a) =Y.[dk + a4'^) Ekk + XY.^>'+I,k, 

k=l k=l 

are equivalent to the following system of differential equations: 

.(1) .(1) 



4 = (l + a4.) 4'-Ci , 



c^P = ci^ (dk+j -dk + ac\}h - ac^^^ 



-k ~ '-k 
■ irn) (m) 



k+j 
(1) 



fe-1 



+ C 



c^ ' ( dk+m - dk + ac 



.(1) 



k+m ^^k-1 I ' 



,(J+1) _ SJ+^)\ 
k ^k-l ) ■> 

1 < j < m — 1. 



Proof - an elementary matrix calculation. 
Notice that 

Bi = 7T+ {{LU-^ - I) /a) , B2 = TT+ {{U-^L - /) /a) 

So, the above Lax representation is indeed of a standard type. 



(15.6) 



(15.7) 



(15. 



(15.9) 



(15.10) 
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One should compare the equations of motion ( [15.9 ) of RTL^^j^(a) with the "first" flow 



of the RTL^_j_^(a) hierarchy (given for m = 2 at the end of Sect. |12D . 

Returning to the relation between the hierarchies RTLj^^^(a) and RTL^^_^(— a), we 
see that it is established via the change of variables 

B(a): 7^Ti^;\(6,a«,...,aM)^7^Tj.:\(d,c«,...,cM 



given by the formulas ( 15. 3| ). Clearly, this transformation is invertible, and its inverse is 
given by 



B-'{a) = B(-a) : 6^ = ^^^, a)^> = ^ . (15.11) 

n (1 + adk+i) 

i=0 



1 + a4 ' '' A 



This transformation relates, as it follows from Theorem 14.2 , the whole hierarchies. In 
particular, this holds for the simplest flows: the change of variables given by the formulas 
(|T5]|) brings the system RTLj;'|^(-a) ( |1LTiD into the system RTlJ^_|^(q) (|15]|), and 



vice versa. 



16 Quadratic r-matrix structure for RTL^^j^(q:) 

Consider the subset 

RT^+i = (^ e go) X I / e g j 

oi g^g, consisting of pairs (L, U^^) of matrices (|15.1| ), (|15.2| ). We show now that this sub- 
set is a Poisson submanifold, ii g0g is equipped with a certain quadratic r-matrix bracket. 



Actually, this bracket was used already in [1^, 16 1 to give an r-matrix interpretation of 



the usual (two-field) relativistic Toda lattice and of the Volterra lattice. 

Theorem 16.1 a) Supply the algebra g ®g with the bracket aPB2(Ai, A2, S) defined by 
the operators 

s -s* ) ' ^* = ( s* -s ) ' ^^^-^^ 



where the operators Ai, A2, S, S* are as in (6J.). Then the set WTm+i is a Poisson 
submanifold in g ® g. 

b) Consider the monodromy maps Mi^2 '■ g ® g ^^ g 

Mi:{L,U-^)^Ti=LU-\ M2 : {L,U-^) ^ T2 = U'^L. (16.3) 

Both maps Mi^2 o,i"c Poisson, if the target space g is equipped with the Poisson bracket 
aPB2(Ai,A2,S). 
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c) Let (p he an Ad-invariant function on g. Then the Hamiltonian equations of motion 
on g (^ g with the Hamilton function if o Mi^2 fnay be presented in the form of the "Lax 
triads " 

L = LC2- CiL, il = UC2 - CiU, (16.4) 

where 

Ci = ^R{dip{LU-^)), C2 = ^R{dip{U-^L)). (16.5) 

Proof. As usual, we have to prove that an arbitrary Hamiltonian vector field on 

(ff®ff,PB2(Ai,A2,S)) 

is tangent to RTm+i- It is enough to consider separately vector fields of Hamilton func- 
tions (f){L) and ^(U), depending only on the first, resp. second, factor in the tensor 
product, because an arbitrary Hamiltonian vector field is a linear combination of such 
ones. Further, the calculations are simplified, if we regard the values of these vector fields, 
parametrizing the second factors via U instead of U^^. With these conventions, the value 
of the vector field of the Hamilton function (piL) in the point oi g ® g parametrized as 
(L, [/), is given by 

2a-^{cp{L),L}2a = L ■ Ai(V0L) - A2{LV^) -L + L- S(LV(/)) - S*(V0L) • L, (16.6) 

2a'^{(p{L),U}2a = U-S*{V(pL)-S{LV(p)-U + U-S{LV(l))-S*{V(l)L)-U. (16.7) 

Similarly, for the vector field of the Hamilton function f{U) in the point oi g g 
parametrized as (L, U), is given by 

2a-^{ip{U),L}2a = -L ■ S{VipU) + S*{U\/ip) -L + L- S{UVip) - S*{VipU) ■ L, (16.8) 
2a-^{ip{U),U}2a = -U-A2{VipU) + Ai{UVip)-U + U-S{U\/ip)-S*{VipU)-U.{16.9) 



Consider the first vector field {(/)(L), •}2a. Its L-component ( 16. 6| ) lies in gig: this is 
a particular case of Theorem |9.1| . The value of its [/-component ( |16.7| ) lies, obviously, in 

m 

Y^ g_j, because the range of the operators S, S* is Qq. Moreover, it is easy to see that 

i=o 

the giQ-component of ( |16.7| ) vanishes. Henceforth, the vector fields {4>{L), ■}2a are tangent 

to WTm+l- 

Let us turn now to the second vector field {(p{U),-}2a- Its L-component, ( [16. 8 ), may 
be rewritten as 

2a-^ip{U),L}2a = LD + DL + W{D)L - LW{D) 
-LD' -D'L-W{D')L + LW{D'), 

where 

D = PoiUVip), D' = Po(V<^[/). 
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The vanishing of the gfj^-coinponent of this expression follows from ( |9.7[) . Finally, for the 
[/-component ( 16.9 ) of the vector field under consideration, we have the following two 
equivalent expressions: 

2a'^{if, U]2a = -2C/P>o(VvjC/) + 2P>o(C/Vc^)C/ 

-C/Po(V(/7C/) + Po(C/Vc/?)C/ + C/Po([/V99) - Po(V(/jC/)[/ (16.11) 

+C/W(Vv?C/) + W{UVip)U - UW{UVip) - W{VipU)U, 



and 



2a-i{^, U}2a = 2t/P<o(V^C/) - 2P^o{UVip)U 

+UPo{VipU) - Po{UVip)U + UPoiUVif) - Po(V99[/)?7 (16.12) 

+UW{VipU) + W{UVip)U - UW{UVip) - W{VipU)U. 



(as compared with the formulas (|9.2| ) and ( |9.3D , respectively, in our present formulas the 

overall signs in the first lines and in the first two terms on the second lines are opposite). 

-1 
The first expression above assures that for t/ E /© gi^ this vector belongs to Qj. 

j=—m i>— "i 

The second expression yields that this vector belongs also to Qj; moreover, since the Qq 

j<o 
component of U is equal to I, it is easy to see that the g^ component of (16.12) vanishes, 

-1 
so that actually it belongs to Qj. Hence, it belongs to Qj. This finishes the 

J<-1 j=-m 



proof of the first statement of the Theorem. The parts b) and c) follow from Theorem 5.1 
with n = 2. I 

Theorem 16.2 The coordinate representation of the bracket induced by aPB2(Ai, A2, S) 
on the submanifold RT^+i, is given by the following formulas: (in ( 16. 13^ ), ({16.16 ) below 
we assume that i < j): 

|4,4'H =-cl^\l + adk), |4'\4+j| =-cl^\l + adk+j), (16.13) 

{4\4^.},„ = -«4^cg, + aci^-^'^c^-p il<l<^-l), (16.15) 

{4^^ 41.},^ = -"4^'^4|. + «4"''^4^f iJ-^ + l<i<J-l). (16.16) 

Proof. The vanishing of all Poisson brackets {dk,d£}2a is a light consequence of ( [16.6 ). 

Similarly, ( 16.13| ) follows simply by considering the Qq component of ( [L6.10| ). The longest 

(7) 

is the calculation of Poisson brackets among the c^ variables. To this end we use the 

defining formula, following from ( |16.11| ): 

2a-^{ip{U),ij{U)}2a = -2(P>o(V(^[/), VipU) + 2{P>o{UVip), UV^j) 
-(Po{VipU),Vi>U) + (Po(t/V<^), C/VV) + (Po(f/Vv9), VVf/) - (Po(V<^C/), UVtP) 
+{W{V^pU),V'tpU) + {W{UV^),UV^) - {W{UVip),W^pU) - {W{VipU),UVip). 

(16.17) 



r-Matrices for Relativistic Deformations of Integrable Systems 443 



Set in this formula 




^{U) = -a4\ 


m) = -«cg„ 


so that 




Vip = X'Ek+^,i, 


S/ijj = X^Ek+e+j,k+i- 


We find: 




VipU = \'Ek+^,k - 


m 




/3=1 



(16.18) 



C/V99 = yEk+^,k - a ^ \'-Pcfl._pEk+i-f3,k. (16.19) 



/3=1 



We consider first the contribution to the Poisson bracket a ^{^,il^}2a = ol< cj^ jCj^^l^ [ 
from the first hne in (|16.17|): 



i— 1 m 



X^Ek+i^k — Oi 2_^ ^^ c^ Ek+i^k+p, a 2_^ X^ '^Cj^j^^Ek+£+j,k+i+j 

\ 13=1 7=i+l / 



i-1 



-iX'-Ek+i^k-a^y ^c,^^._^Ek+i-f3,k,a ^ A^ '^CiJ^g^j_^Ek+i+j-y,k+i) ■ 

\ /3=1 7=i+l / 

Calculating these scalar products, we find: 

m 

= a 22 [ck+i^e+j^iZ+i - '^k+e+j-'y^o,i+j-y^i,eJ 

7=j+l 

i— 1 m 

13=1 7=j+l 

~ "^fc+^ "^-3 "'^fc "^,1 " ^ke. ^k '-k+e ' 
where 

^M =xK-i + l;*-i- '^)xdi-m,i-j - 1) -xei'i^^m- j)xiil,i- 1) ,^^ ^^, 

(lD.2Uj 

= X£(-min(j - l,m-i),i- j - 1) - x<? (l,min(i - l,m- j)) . 

Assuming, for the sake of definiteness, that i < j, we see that the intervals of the two 
characteristic functions in the last line do not intersect. So, we found the contributions to 
the Poisson bracket a < ci, ,cf!^^ \ described by the first term on the right-hand side of 
the formula ( 16.14 ), and the second terms on the right-hand sides of ( 16.15| ), ( |16.1(j ). 



Finally, calculating the contribution of the remaining part of ( |16.17D , we get: 
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where 



Mj) 



'M 



+Wk+i,k+e+j + Wk,k+e - Wk,k+e+j — Wk+i^k+e)- 
This coefficient differs from ( p. 21 ) by —6i^£^j + (5o/, so that the result (for i < j) 
written down immediately: 

^kP = Xi{-j, -J + ^ - 1) - X£(l, i)- 
This describes the first terms on the right-hand sides of ( 16.15| ), ( |16.16 ), which finishes 
the proof of the part a) of Theorem 16.1 and the proof of Theorem 16.2 . The parts b) 
and c) of Theorem 16.1| follow from Theorem fi.l\ with m = 2. I 

Corollary. The system RTLj^_j_-^(a) UlS.dj ) is Hamiltonian with respect to the bracket 



(16.22) 
may be 

(16.23) 



(16.1S)-(16.16), with the Hamilton function 

N N 



a 



I))o- 



(16.24) 



fc=i 



fc=i 



(+) 



Unfortunately, the RTL^_j|.j^(a) hierarchy seems not to have an invariant linear Poisson 
bracket. Nevertheless, it is bi-Hamiltonian, and the most direct way to find its second 
Hamiltonian formulation is to push forward under B(a) the known invariant Poisson 

brackets of the RTL^^]^(— a) hierarchy. It can be verified that the bracket {•, •}2a found 
in Theorem 16.2 is recovered by applying this construction to the bracket {•, ■}i — a{-, ■}2 



(-) 



of the hierarchy RTL^_j^]^(— a). A compatible bracket appears, if we start just with {•, ■}i. 
We formulate this as a separate statement. 

Theorem 16.3 The hierarchy RTL^_j_-^(a) has an invariant cubic Poisson bracket {•, -y^a) 
which can be defined as the push-forward under the change of variables B(a) ( fJ5. j| j of the 
invariant bracket {•,•}! of the RTL^_j_-^(— a) hierarchy. The brackets {•,-}2a and {^-jsct 
are compatible. 

We shall not try to write down the general formulas for this pushed bracket, as they 
turn out to be complicated enough. However, in the next Section we give the result of 
calculations for this bracket in the particular case m = 2. 



17 Example: RTL^+^(q:), 

the three-field analog of the relativistic Toda lattice 

In order to provide the reader with an illustration of the bi-Hamiltonian structure of the 



system RTL^_j^^(a), we give here the corresponding formulas for the case m = 2, 
is the next in complexity case after the usual relativistic Toda lattice RTL4_(a) 

equations of motion of the system RTLg (a) read: 



which 
, The 



4 = (1 + adk) 






a(1) 



a(2) 



4 ( 4+1 -dk + ac 



(1) 
fc+i 



.(2) 



dk+2 



dk + "4+2 



ac 



ac 



(1) 
fe-i 

,(1) 
fc-i 



+ 



,(2) _ „(2) 
'k ^fc-1 



(17.1) 
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This is, of course, an 0(a)-perturbat 


ion of the system (10.1). The quadratic Poisson 


structure of Theorem 16.1 is charachterized by the following nonvanishing brackets: 


1- J 2a 


{4'\4+i}^^ = -4'^(l + a4+i), 


{4,4'^}^^ = -cf(l + a4), 


{4'\4+2}^^ = -4'^(l + a4+2), 




r j2) (1) 1 _ __(2) (1) 


r J2) J2) 1 _ _^J2) (2) 


1^(2) ^(2) 1 _ _„c(2) (2) 





This quadratic bracket is an 0(a)-perturbation of the linear invariant bracket (10.2) of 
the nonrelativistic system TL3. The cubic invariant Poisson bracket (see Proposition 16.3 ) 
is characterized by the fohowing nonvanishing brackets: 

{4, 4+i}3a = -"4^'' (1 + «^fc) (1 + adk+i) , 



dk, 4'^}3^ = -4'^ (1 + «4) (1 + adk + a' 4'^) ' 

4 , 4+1 1 = -4 (1 + «4+i) (1 + a4+i + "^4 ) ' 

4, 4+1)3^ = -« (4^^ + «4^^4+i) (1 + «^fc) ' 



'-k ' 



4+2I = -a (^4 + «4 4+1) (1 + «4+2) , 



^fc'4 } = -c[' (l + adk) \l + adk + a^cl' 



k ' 



^fc+if, =-« 4 (i + «^fe+i)(4 -Cfc+i), 



3a 



4, 4+1)3^ = -«^4^^4+i (1 + «4) 



'-k ' 






4+2 [ = -4 (1 + «4+2) ( 1 + a4+2 + a^4+i 
4+3I = -a^4^^4+2 (1 + "4+3) , 



Ji) Ji) \ _ {^,^^ ^ Z'^^^) „ (1) (1) 



(17.3) 



\^fc '^fc+2/ 



2 (2) (1) 2 (1) (2) 3 (1) (1) (1) 

3^ ~ " ^fc ^fc+2 " '-k ''fc+1 " ^fc ^fc+l^fc+2' 
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4^^ 4^i}3„ = — r^4^i (1 + «4.o - a^4^) (4^^ + -c^^'c^l^) , 

4^^ 41}3„ = -2«4^^4^i (1 + «4.o - a^4?, (4^^ + a4^) (4^^ + 4^i)) ' 

4^^ 4^-J.}3„ = -2-f 4^i (1 + «4..) - a^4^) (cS, + acg, (4^), + ^ 

'-fc '^A:+2J3^~ " ''fc+2 l^^fc T^ "^A: ^k+1 



i^\ 4+1} = -"4^^4+1 (2 + adk+1 + a4+2 + a^4^^ + «^4 



f ^ 4+2)3^ = -«4'^4+2 (2 + 2a4+2 + «'4+i + «'4+2 

.(2) J2) "1 _ _^3 J2) (1) (2) 
'fc '''fc+3/3^~ " ''fc ''fc+2''fc+3' 

This cubic bracket is again an 0(a)-perturbation of the linear invariant bracket ( |10.2| ) of 
the nonrelativistic system TL3. However, the Unear combination 

a~^ {{■,-}3a - {•,-}2a) 

of the both compatible Poisson brackets above leads, in the limit a ^ 0, to the quadratic 
invariant bracket (10.3) of TL3. 



18 Conclusions 

The results of this paper confirm once more that the r-matrix theory, with its various 
generalizations, provides very convenient means for studying the Hamiltonian aspects of 
integrable lattice systems. The relativistic lattice KP, which seemed for a long time to 
lie outside of the applicability area of the r-matrix theory, actually fits very nicely into 
this framework. Moreover, it allowed us to find for the first time the quadratic invariant 
bracket of the relativistic lattice KP, and to establish an amasing fact: that this attribute 
remains undeformed by the relativistic deformation. By the way, this gives also an answer 
to a problem left open in [O], concerning the Hamiltonian nature of the "relativistic 
Bogoyavlensky lattices". Namely, these relativistic systems are Hamiltonian with respect 
to the invariant quadratic brackets of their nonrelativistic counterparts. A detailed account 
of this latter result will be given elsewhere. 
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